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ABSTRACT 

A  concept  for  a  terahertz  laser  in  vapor -phase-grown  homoepitaxial  GaAs  with  spatially  periodic  doping  profde  was 
theoretically  explored.  Monte  Carlo  simulation  of  hole  transport  in  multilayer  delta-doped  p-GaAs/GaAs  structures  in 
crossed  electric  and  magnetic  fields  was  performed  to  investigate  possibilities  of  the  terahertz  amplification  on  inter¬ 
valence-hand  light-to-heavy  hole  transitions.  The  results  are  compared  to  those  calculated  for  uniformly  doped  bulk  p- 
GaAs  and  recently  proposed  p-Ge/Ge  structures.  The  improvement  in  the  gain  for  delta-doped  p-GaAs  structures  is 
about  ~  2  -  3  times  over  bulk  p-GaAs.  Terahertz  laser  generation  in  the  considered  GaAs  device  concept  appears 
feasible,  as  is  growth  of  structures  with  active  thicknesses  sufficient  to  support  quasioptical  cavity  solutions  at  100  pm 
vacuum  wavelengths.  Potential  applications  for  the  considered  laser  device  include  sensing  of  chem/bio  agents  and 
explosives,  biomedical  imaging,  non-destructive  testing,  and  communications. 

Keywords:  terahertz,  far  infrared,  laser,  germanium,  silicon,  gallium  arsenide,  epitaxy. 


1.  INTRODUCTION 

The  THz  region  of  the  electromagnetic  spectrum,  loosely  defined  here  to  span  the  range  0.3  to  20  THz,  is  of  particular 
interest  to  various  department  of  defence  agencies.  Imaging  through  dust,  or  so  called  “brown-ouf  ’  conditions,  as  well 
as  through  fog,  are  of  particular  interest  to  both  the  Army  and  the  Air  Force  for  various  tactical  reasons.  In  such  cases 
there  is  always  a  tradeoff  between  radiation  loss  and  resolution.  Longer,  mm-waves  have  low  loss  but  poor  resolving 
power,  while  shorter,  IR  waves  produce  highly  resolved  images  but  with  extremely  poor  penetration  depth.  This 
problem  of  loss  versus  resolution  motivates  investigations  into  the  lower  end  of  the  above  mentioned  spectral  range  with 
the  added  constraint  that  any  such  system  must  also  be  portable  in  order  to  be  of  use  in  the  field.  The  mid-THz  range 
would  show  increased  resolution  when  used  in  an  imaging  system  but  because  of  the  increased  loss  is  more  suitable  for 
use  in  well  controlled  environments  such  as  personnel/package,  and  non-destructive  industrial  screening,  as  well  as  non- 
evasive  medical  imaging.  Beyond  imaging  systems,  the  Air  Force  is  interested  in  the  mid-  to  upper-end  of  the  THz 
spectrum  for  use  in  advanced  threat  warning  systems,  where  biological  and  chemical  agents  have  rich  spectra,  as  well  as 
remote  earth  sensing,  and  flame/plasma  characterization. 

Based  on  the  above  needs,  a  great  interest  has  arisen  in  developing  compact,  portable,  THz  sources.  At  present,  because 
the  THz  spectrum  constitutes  a  gap  between  well-developed  electronic  and  photonic  source  technology,  both  electronic 
and  photonic  approaches  have  been  extended  to  produce  THz  radiation.  Frequency  multiplied  electronic  sources  have 
demonstrated  3  mW  of  1.9  THz  radiation  with  just  0.2%  conversion  efficiency,*  while  quantum  cascade  lasers  have 
demonstrated  4  mW  of  continuous-wave  radiation  at  4.4  THz  with  33%  efficiency  at  an  operating  temperature^  of  48  K. 
Non-linear  optical  techniques,  such  as  optical  recification  (OR)  have  also  be  used  for  THz  generation.  Although 
chromatic  dispersion  effects  in  a  highly  non-linear  crystals  such  as  GaAs,  have  historically  limited  the  use  of  this 
material  in  OR  based  systems,  quasi-phase-matching  has  recently  been  used  to  counter  dispersion  effects  and  produce 
THz  radiation.^  In  Ref  3,  quasi-phase-matching  was  achieved  through  the  growth  of  thick  GaAs  layers  via  hydride 
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vapor  phase  epitaxy  (HVPE)  on  orientation-patterned  substrates/  The  development  of  the  HVPE  teehnique,  eapable  of 
produee  epilayers  of  GaAs  100s  of  mierons  thiek,  along  with  the  potential  to  delta-dope  this  material  p-type,  motivates 
the  investigation  of  this  material  for  use  as  a  THz  gain  medium  utilizing  inter-valenee-band  (IVB)  effeets  in  a  fashion 
similar  to  what  has  reeently  been  proposed  for  p-Ge/Ge  struetures/’^ 

As  previously  stated,  GaAs  is  a  potential  aetive  medium  for  inter-valenee-band  (IVB)  hot-hole  THz  lasers/'^  Monte 
Carlo  simulations  have  shown  that  bulk  uniformily-doped  GaAs  should  have  aeeeptable  performanee.  However,  due 
to  doubled  (see  below)  impurity  seattering  rate  (a  eonsequenee  of  lower  dieleetrie  eonstant)  in  GaAs  relative  to  Ge,  THz 
gain  in  p-GaAs  deviees  is  expeeted  to  be  lower  than  in  the  well-established  p-Ge  lasers.  Compared  with  Ge,  GaAs  has 
additional  faetors  of  polar  optieal  phonon  and  aeoustie  piezoeleetrie  phonon  seattering  proeesses.  The  latter  is  usually 
unimportant,  but  the  polar  optical  phonon  scattering  has  a  significant  role.  Monte  Carlo  simulation  results  for  GaAs 
illustrate  these  features. 

For  the  IVB  mechanism,  the  inverted  population  grows  at  certain  ratios  of  applied,  crossed,  electric-  and  magnetic- 
fields,  when  light  holes  are  accumulated  on  closed  trajectories  below  the  optical  phonon  energy,  while  heavy  holes 
undergo  rapid  optical  phonon  scattering.  Fig.  1  is  a  schematic  quasiclassical  explanation  of  the  IVB  laser  mechanism 
that  is  valid  at  low  applied  fields  where  Landau-level  spacing  is  much  less  than  carrier  kinetic  energy.  The  emission  of 
polar  (circle)  and  non-polar  (starburst)  optical  phonons  is  indicated.  Polar  optical-phonon  scattering  in  III-V 
compounds’*  provides  a  “hard  ceiling”  for  hot  holes.  Most  of  the  scattering  events  occur  right  beyond  the  optical 
phonon  threshold,  which  is  35  meV  for  GaAs.’^  In  contrast,  in  non-polar  semiconductors  (Ge,  Si)  the  optical-phonon 
emission  rate  grows  slowly  with  energy  from  the  threshold  value.  Light-hole  life  time,  which  is  responsible  for  the 
inversion  population  and  for  the  gain,  is  determined  by  acoustic  phonon  scattering,  ionized  impurity  scattering,  and 
carrier-carrier  interaction.  This  limits  the  operation  of  the  IVB  hot-hole  laser,  which  has  been  experimentally  realized 
only  for  p-Ge,  to  low  temperatures  (T  <  20  K)  and  low  carrier  concentration  (p  ~  lO’"*  cm'^).’"^  The  IVB  amplification 
mechanism  illustrated  in  Fig.  1  provides  wide  tunability  of  p-type  hot-hole  devices  in  the  terahertz  frequency  range  (1.5 
-  4.2  THz  for  p-Ge). 


Fig.  1.  Schematic  of  quasi-classical  explanation  for  intersubband  hot-hole  laser. 


We  recently  developed  a  terahertz  laser  concept  based  on  inter-valenee-band  transitions  of  holes  with  transport  in 
crossed  electric  E  and  magnetic  B  fields  in  a  planar  periodically  doped  p-Ge/Ge  structure.^'^  The  design,  shown 
schematically  in  Fig.  2,  achieves  spatial  separation  of  hole  accumulation  regions  from  the  doped  layers,  which  reduces 
ionized-impurity  scattering  and  carrier-carrier  scattering  for  the  majority  of  light  holes,  allowing  significant  increase  of 
total  carrier  concentration  without  affecting  light  hole  life  time.  The  resulting  increase  in  gain  over  the  bulk  p-Ge  laser 
promises  to  raise  maximum  operation  temperatures  to  77  K.  At  the  same  time  the  proposed  laser  retains  the 
intersubband  mechanism  with  its  wide  tuning  range  1  -4  THz.  Moreover,  this  crystalline-Ge  device  can  be  grown  by 
chemical  vapor  deposition  (CVD),  which  allows  active  thicknesses  comparable  to  the  THz  wavelength,  thus  allowing 
low-loss  quasi-optical  cavity  solutions.  As  soon  as  we  have  developed  a  THz  laser  concept  based  on  CVD  epitaxial  Ge 
devices,  it  becomes  very  interesting  to  consider  the  possibility  of  realizing  a  similar  structure  based  on  GaAs,  because  of 
the  highly  sophisticated  technology  available  for  epitaxial  growth  of  this  material. 
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Fig.  2.  FLomoepitaxial  multilayer  p-type  semiconductor  structure  TFIz  laser  concept  for  germanium.  The  layer  period  is  chosen  to  be 
200  -500  nm,  which  is  larger  than  the  light  hole  cyclotron  orbit,  but  smaller  than  that  for  heavy  holes.  The  doped  layers  are  ~10%  as 
thick  as  the  period.  Light  holes  are  shown  by  Monte  Carlo  simulation  to  accumulate  in  the  undoped  layers  where  their  life  time  is 
enhanced  by  reduced  impurity  scattering.  This  scheme  allows  much  higher  average  carrier  concentrations  than  for  uniformity  doped 
bulk,  giving  a  corresponding  increase  in  gain.  The  stack  thickness  achievable  for  germanium  by  chemical  vapor  deposition  is 
estimated  to  be  100  |am.  In  the  case  of  a  GaAs  structure,  the  total  stack  thickness  that  can  been  achieved  by  vapor  phase  epitaxy 
would  be  up  to  800  microns,  giving  a  low-loss  quasioptical  cavity  solution. 

2,  THEORETICAL  METHODS 

Hole  dynamics,  hole  distribution  functions,  and  the  gain  on  direct  optical  light-to-heavy  hole  transitions  in  p-GaAs  are 
calculated  by  the  Monte-Carlo  simulation  method  using  classical  motion  equations  and  hole  scattering  probabilities.’^  ’"^ 
Two  valence  subbands  (light  and  heavy  holes)  with  isotropic  and  parabolic  dispersion  laws  are  considered.  The 
isotropic  approach  is  justified  by  the  relatively  small  warping  of  the  GaAs  valence  band.  We  can  neglect  quantum 
confinement  effects  because  the  considered  structure  has  no  heteroboundaries,  only  delta-doped  layers.  We  also  neglect 
Landau  quantization,  which  is  a  good  approximation  for  magnetic  fields  of  3  -  4  T  and  below.  Time  or  ensemble 
averaged  momentum  and  position  yield  the  hole  distribution  functions  fi,h(k,r)  (subband  (l,h),  wavevector  k,  and 
coordinate  r).  The  distribution  functions  were  considered  uniform  in  the  horizontal  planes  according  to  the  geometry  of 
the  problem.  The  standard  Rees  rejection  technique  chooses  among  scattering  processes.’^  The  rate  of  each  scattering 
process  is  given  by  a  temperature-dependent  analytic  expression.  Polar  and  non-polar  optical  phonon  scattering  is 
treated  in  a  deformation  potential  approximation.”  Acoustic  phonon  scattering  is  simplified  according  to  Ref  15. 
Inelasticity  for  acoustic  phonon  scattering  is  included.’^  Piezoelectric  scattering  is  treated  according  to  Ref  17.  Brooks- 
Herring  model”  with  inverse  Debye  screening  length  and  Yukawa  potential  was  used  for  ionized  impurity  scattering. 
Hole-hole  scattering  was  calculated  according  to  Ref  18.  Iteration  determines  the  self-consistent  solntion  of  the 
Poisson  equation  and  thereof  the  spatial  carrier  distribution  and  potential  profile.  The  small  signal  gain  is  calculated  as 
the  difference  between  the  gain  on  direct  inter-valence-band  (light  to  heavy  hole)  transitions  and  free  carrier  absorption 
assisted  by  phonons  and  ionized  impurities. Up  to  liquid  nitrogen  temperatures  and  impurity  concentrations  of  ~  lO’^ 
-  lO’^  cm”,  the  main  contribution  to  free  carrier  absorption  in  GaAs  comes  from  polar  optical  phonon  emission  by  heavy 
holes.  Lattice  absorption^’  was  not  included  in  the  calculations. 

Strong  inter-valence-band  (heavy-to-light)  polar-optical  phonon  scattering  in  p-GaAs  populates  the  light-hole  band.  In 
contrast  to  non-polar  optical  phonon  scattering,  most  of  these  scattering  events  occur  right  beyond  the  sharp  energy 
threshold  for  polar-phonon  scattering  (see  Fig.  1).  As  a  result,  after  emission  of  an  optical  phonon  light  holes  are  created 
with  very  small  kinetic  energy.  These  slow  holes  experience  very  strong  ionized  impurity  and  hole-hole  scattering, 
which  is  enhanced  in  GaAs  compared  to  Ge  because  of  the  smaller  relative  permittivity  {e^  =  12.9  for  GaAs  vs  16.0  for 
Ge).  These  holes  are  quickly  accelerated  by  the  applied  electric  field,  but  they  return  to  the  low-energy  range  on  each 
loop  of  cyclotron  trajectory  in  crossed  fields.  The  ionized  impurity  scattering  is  the  dominant  light-hole-band 
depopulating  mechanism.  Similarly,  ionized  impurity  scattering  of  heavy  holes  is  also  stronger  in  GaAs  compared  to 
Ge.  Acoustic  deformation  potential  scattering  rates  in  GaAs  are  similar  to  rates  in  Ge.  Piezoelectric  scattering  is  weak 
and  can  be  neglected.  In  order  to  maintain  an  inverted  distribution  of  holes  and  to  reduce  the  effect  of  ionized  impurity 
scattering  for  light  holes  in  GaAs,  the  magnitude  of  the  applied  fields  must  be  higher  than  typical  for  the  p-Ge  laser  (E  = 
1  -  2  kV/cm,  B  =  1  -  2  T). 
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3,  CALCULATIONS  FOR  UNIFORMLY  DOPED  MATERIALS. 


For  direct  comparison  of  the  scattering  effects  in  germanium  and  GaAs,  the  distribution  functions  averaged  over 
directions  are  plotted  in  Fig.  3  for  both  materials  in  strong  crossed  electric  (E  =  4  kV/cm)  and  magnetic  (B  =  3  T)  fields 
at  T  =  20  K  for  carrier  concentration  p  =  2.5  x  lO’"^  cm'^.  Note  in  Fig.  3  the  “hard  ceiling”  for  heavy  hole  energies  in  the 
case  of  GaAs  that  occurs  at  the  35  meV  optical  phonon  energy.  This  effect  causes  the  hole  distribution  for  GaAs  to  be 
shifted  into  low-energy  regions  compared  with  Ge. 


Fig.  3.  Calculated  distribution  functions  of  light  and  heavy  holes  in  Ge  and  GaAs.  E  =  4  kV/cm,  B  =  3  T,  T  =  20  K,  and  p  =  2.5x10'"* 
cm'^. 


Calculated  gain  for  amplification  on  direct  inter-valence-band  transitions,  and  free-carrier  absorption,  for  bulk  p-Ge  and 
p-GaAs  are  presented  in  Figs.  4  and  5,  respectively.  The  simulation  parameters  for  Ge  (GaAs)  were  E  =  1.5  kV/cm,  B  = 
1 .2  T,  T  =  lOK,  p  =  lO’"*  cm'^  (E  =  4  kV/cm,  B  =  3.5  T,  T  =  lOK,  p  =  2\  lO’"*  cm'^).  The  total  gain  is  smaller  for  GaAs 
than  for  Ge,  despite  the  much  higher  applied  fields  and  twice  higher  hole  concentration.  Fields  need  to  be  increased  over 
those  for  Ge,  in  principle,  to  achieve  any  positive  gain  in  GaAs.  One  reasons  for  this  is  the  almost  twice  smaller 
oscillator  strength  for  GaAs.^^  A  second  reason  is  stronger  free  carrier  absorption,  which  arises  because  phonon 
scattering  rates  are  higher,  and  this  factor  is  especially  unfavorable  considering  that  the  inversion  between  light-  and 
heavy-hole  bands  is  shifted  to  lower  energy  in  GaAs. 


Fig.  4.  Calculated  inter-valence-band  gain  and  free  carrier  absorption  in  bulk  p-Ge  with  carrier  concentration  lO'"*  cm‘^.  E  =  1 .5 
kV/cm,  B  =  1 .2  T,  T  =  lOK.  Total  gain  (smooth  curve)  is  the  difference  between  two  curves  with  symbols. 
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Wavenumber  (cm'^) 

Fig.  5.  Calculated  inter-valence-band  gain  and  free  carrier  absorption  in  bulk  p-GaAs  with  carrier  concentration  2x10'^'  cm‘^.  E  =  4 
kV/cm,  B  =  3.5  T,  and  T  =  10  K.  The  total  gain  (smooth  curve)  is  the  difference  between  two  curves  with  symbols. 

The  maximum  value  of  the  ealeulated  gain  in  bulk  p-GaAs  is  presented  in  the  eontour  plot  Fig.  6  as  a  funetion  of  the 
applied  fields.  The  interval  between  eontour  plots  is  0.025  em  ’.  Flere  and  below,  the  contours  are  spline-interpolated. 
The  waves  in  the  contours  are  caused  by  the  interpolation  to  the  limited  set  of  data  points,  which  occur  at  every  grid 
intersection.  As  has  been  shown  previously, the  optimal  ratio  of  fields  and  their  magnitudes  must  be  different  for  the 
two  materials,  and  the  gain  in  GaAs  can  be  increased  by  applying  higher  electric  (~  8  kV/cm)  and  magnetic  (~  7  T) 
fields.  Flowever,  the  validity  of  classical  approach  for  Monte  Carlo  simulation  under  these  extreme  conditions  becomes 
questionable.  Performance  of  p-GaAs  hot-hole  lasers  would  be  improved  by  the  periodic  doping  scheme,  suggested  in 
the  introduction,  by  which  impurity  scattering  is  eliminated  from  the  active  region.  The  same  device  concept  with 
vertical  transport,  which  has  been  demonstrated  theoretically  for  Ge,^  can  be  applied  directly  to  GaAs.  This  will  be 
discussed  next. 


Fig.  6.  Gain  for  uniformly  doped  p-GaAs  as  a  function  of  applied  fields.  Simulation  parameters  are  p  =  2x10*"^  cm'^,  T  =  10  K. 
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4,  CALCULATIONS  FOR  DELTA  DOPED  STRUCTURE 


A  delta  doped  p-GaAs/GaAs  structure,  with  doping  period  d=  200  nm  (20%  doped,  80%  undoped)  and  average  carrier 
concentration  2x10*"^  cm'^,  was  used  for  preliminary  gain  estimations.  The  electric  and  magnetic  fields  are  4  kV/cm  and 
3.5  T,  respectively,  and  the  lattice  temperature  is  10  K.  Calculated  spatial  distribution  of  light  and  heavy  holes  across 
two  structure  periods  is  presented  in  Fig.  7.  Fleavy  holes  have  a  slightly  larger  concentration  near  the  doped  regions. 
Most  of  the  light  holes  are  concentrated  in  the  undoped  region  between  two  doped  layers.  Because  of  the  much  larger 
density  of  states  of  the  heavy  hole  band,  the  calculated  distribution  in  the  undoped  regions  constitutes  a  population 
inversion. 


Fig.  7.  Distribution  of  light  and  heavy  hole  concentrations  across  the  structure.  E  =  4  kV/cm,  B  =  3.5  T,  <7  =  200  nm,/)  =  2x10*"*  cm‘^. 

The  resulting  spatial-spectral  gain  distribution  is  shown  in  Fig.  8.  Strong  inversion  is  observed  in  the  regions  of  undoped 
GaAs.  The  doped  layers  are  absorbing. 


3.5T,  T  = 


The  spatially  averaged  net  gain  in  the  structure  is  presented  in  Fig.  9.  Compared  to  Fig.  5  for  bulk  p-GaAs  the  gain  is 
approximately  twice  higher.  The  gain  spectrum  is  ~  50  -  150  cm  *  wide,  but  this  will  be  modified  by  the  lattice 
absorption  (not  included  in  the  calculations).  Note  that  all  simulation  parameters  are  the  same  as  in  Fig.  5  for  bulk  GaAs 
(fields,  average  concentration,  temperature)  except  the  doping  profile,  so  that  the  observed  improvement  is  purely  an 
effect  of  the  selective  doping.  As  will  be  shown  below,  an  additional  improvement  of  the  gain  by  increase  in  the  carrier 
concentration,  as  was  observed  for  p-Ge/Ge  structures,^’^  is  not  achievable  for  GaAs. 
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Fig.  8.  Spatial-spectral  gain  distribution  across  two  structure  periods  (Nav=  2  x  lO'^'cm’^,  d=  200  nm,  E  =  4  kV/cm,  B  ^ 
10  K)  with  free  carrier  absorption  included. 
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Fig.  9.  Calculated  spatially-averaged  terahertz  gain  and  free  carrier  absorption  in  the  delta-doped  GaAs  structure  with  all  other 
parameters  the  same  as  for  the  uniformly  doped  case  presented  in  Fig.  5.  The  total  gain  (smooth  curve)  is  the  difference  between  two 
curves  with  symbols. 

Next,  we  optimize  the  strueture  by  varying  the  simulation  parameters  to  maximize  the  gain.  The  dependenee  of  peak 
gain  on  the  magnitude  of  the  applied  fields  for  the  strueture  with  fixed  doping  parameters  (p  =  2xl0’"^  em'^,  d=  250  nm, 
10%  doped,  90%  undoped)  is  presented  in  the  eontour  plot  Fig.  10.  The  interval  between  eontours  is  0.025  em  '. 
Comparison  of  Figs.  6  and  10  shows  that  the  aeeeptable  range  of  the  applied  fields  is  wider  for  the  delta-doped  strueture 
than  for  the  uniformly  doped  GaAs  (see  0.15  em  ’  eontours  in  both  figures).  The  gain  aehievable  for  delta  doped 
stmeture  is  about  twiee  higher  than  for  bulk  p-GaAs. 


Magnetic  field  (T) 

Fig.  10.  Gain  for  delta-doped  p-GaAs/GaAs  structure  as  a  function  of  applied  fields.  Simulation  parameters  are  p  =  2x10'"*  cm‘^,  d  = 
250  nm  (10%  doped,  90%  undoped),  T  =  10  K. 
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Next,  we  eonsider  the  effeet  of  varying  the  doping  parameters  (strueture  period  and  average  earrier  eoneentration)  while 
holding  the  applied  fields  fixed.  These  fields  (E  =  4.5  kV/em,  B  =  4  T)  eorrespond  to  the  point  indieated  by  an  asterisk 
in  Fig.  10.  The  relative  thiekness  of  the  p-GaAs  doped  layers  is  fixed  and  equal  to  10%  of  the  period  as  before.  The  gain 
surfaee  is  shown  in  Fig.  11  by  eontours  with  0.025  em  ’  intervals.  In  eontrast  to  p-Ge/Ge  struetures,^  the  optimal 
average  earrier  eoneentration  for  the  delta  doped  p-GaAs/GaAs  strueture  is  eomparable  to  that  for  uniformly  doped  p- 
GaAs’°  and  equal  to  ~2  x  lO’"^  -  4  x  lO’"^  em'^.  The  optimal  range  of  the  strueture  period  is  250  -  500  nm  for  the  ehosen 
fields.  As  ean  be  seen  from  the  Fig.  11,  the  gain  is  very  sensitive  to  the  strueture  period  at  high  doping  eoneentrations. 
This  is  eaused  by  the  distortion  of  eleetric  potential  across  the  structure  period,  which  changes  the  ratio  of  the  local 
electric  and  magnetic  fields  and  hole  lifetime.  At  low  carrier  concentration  the  gain  surface  does  not  show  such  a  strong 
dependence  on  structure  period  as  soon  as  the  average  light  hole  cyclotron  orbit  (which  decreases  with  increase  in 
magnetic  fields)  fits  inside  the  undoped  layer  (see  decrease  in  the  gain  in  the  low  period  part  of  the  Fig.  1 1)  to  provide  a 
long  lifetime  of  the  light  hole  and  therefore  inversion.  The  asterisk  in  Fig.  1 1  indicates  the  structure  parameters  for  the 
calculations  in  Fig.  10. 
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Fig.  1 1 .  Gain  for  p-GaAs/GaAs  structure  as  a  function  of  carrier  concentration  and  structure  period.  Simulation  parameters  E 
kV/cm,  B  =  4T,  T  =  10  K.  The  relative  thickness  of  the  doped  p-GaAs  layers  is  10%  of  the  structure  period. 
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5,  GROWTH  FEASIBILITY 

Thick  epitaxial  GaAs  is  made  possible  by  hydride  vapor  phase  epitaxy  (HVPE),  which  can  achieve  much  faster  growth 
rates  than  the  more  widespread  GaAs  epitaxial  growth  methods  (molecular  beam  epitaxy  or  organometallic  vapor  phase 
epitaxy).  We  have  demonstrated  growth  rates  well  above  100  pm/hr  for  homoepitaxial  GaAs  via  HVPE  in  a  custom- 
built  reactor  at  the  Air  Force  Hanscom  Research  Site.^^  The  system  consists  of  a  horizontal  quartz  tube,  heated  by  a 
three-zone  furnace  and  sealed  to  allow  low  pressure  operation  in  the  range  of  1  to  5  torr  (Fig.  12).  HCl  vapor  passing 
over  a  liquid  Ga  source  reacts  to  form  GaCl,  which  is  transported  to  the  substrate.  Arsenic  is  supplied  in  the  form  of 
arsine  (ASH3),  which  decomposes  on  the  surface.  Both  the  Ga  source  and  the  substrate  are  heated  to  temperatures  in  the 
range  of  650  to  750°C;  the  temperature  gradient  along  the  tube  drives  GaCl  formation  at  the  source  and  GaAs  deposition 
on  the  substrate.  The  growth  rate  can  be  controlled  by  varying  the  vapor  supersaturation  -  i.e.,  increasing  or  decreasing 
the  partial  pressures  of  the  reactants  with  respect  to  their  equilibrium  values  alters  the  tendency  toward  growth  or 
etching  of  GaAs.  This  control  allows  us  to  achieve  very  rapid  epitaxial  growth  rates,  and  we  regularly  achieve 
thicknesses  in  the  range  of  500  to  800  pm  in  single  HVPE  runs.  The  p-type  dopant  will  be  introduced  through  an 
additional  port  and  mixed  with  the  gases  upstream  of  the  wafer.  At  a  growth  rate  of  100  pm/hr,  a  doping  period  of  200 
nm  corresponds  to  8  seconds.  Abrupt  doping  transitions  are  facilitated  in  our  system  by  the  use  of  a  rapid  gas  velocity 
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through  the  reactor  at  low  pressure,  allowing  for  very  short  dopant  residence  times.  Therefore,  the  growth  of  multilayer 
delta-doped  p-GaAs/GaAs  structures  of  up  to  800  pm  total  thickness  is  attainahle. 


Fig.  12.  Schematic  of  low-pressure  hydride  vapor  phase  epitaxial  reactor.  GaAs  growth  rates  of  up  to  150  pm/hr  and 
total  thicknesses  of  800  pm  have  heen  demonstrated. 


6,  SUMMARY 

Gain  calculations  for  a  delta-doped  homoepitaxial  p-GaAs/GaAs  THz  laser  structure  were  performed,  and  the  results 
were  compared  with  recently  published  calculations  for  a  similar  p-Ge/Ge  laser  structure.  Though  the  GaAs  structure 
gives  lower  gain  and  requires  higher  applied  fields  than  found  for  Ge,  GaAs  has  the  advantage  of  the  existence  of  more 
highly  developed  epitaxial  growth  technology.  Indeed,  consideration  of  the  potential  of  a  particular  vapor-phase-epitaxy 
growth  reactor  suggests  that  growth  of  the  considered  device  structure,  with  sufficient  active  thickness  for  low-loss 
quasi-optical  cavity  solutions,  is  feasible.  What  remains  to  be  done  is  experimental  demonstration  of  the  growth  of  the 
device  with  comparison  of  physical,  electrical,  and  optical  characterization  data  to  design  specifications  and  comparison 
with  calculations.  Expectations  are  the  development  of  a  new  type  of  THz  laser  technology  that  can  compete  favorably 
with  MBE  grown  quantum  cascade  lasers. 
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